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Abstract 

We calculate power spectra of cosmological perturbations at high accuracy 
for two classes of inflation models. We classify the models according to 
the behaviour of the Hubble distance during inflation. Our approximation 
works if the Hubble distance can be approximated either to be a constant 
or to grow linearly with cosmic time. Many popular inflationary models can 
be described in this way, e.g., chaotic inflation with a monomial potential, 
power-law inflation and inflation at a maximum. Our scheme of approxi- 
mation does not rely on a slow-roll expansion. Thus we can make accurate 
predictions for some of the models with large slow-roll parameters. 
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1 Introduction 



The most important prediction of cosmological inflation jl] (besides the spa- 
tial flatness of the universe) is the generation of primordial fluctuations of 
matter and space-time [0. These fluctuations seed the formation of large 
scale structure and give rise to anisotropics in the cosmic microwave back- 
ground (CMB). Recent CMB observations || have clearly detected the first 
three acoustic peaks. In order to explain these observations a dominant 
contribution of adiabatic perturbations is necessary, as predicted by cosmo- 
logical inflation. The position of these peaks is consistent with a spatially 
flat universe. The level of accuracy of current experiments is of the order of 
10%, which roughly corresponds to 10% uncertainty for the best determined 
cosmological parameters. Future experiments [[| will increase this accuracy 
to the limit of the cosmic variance (of the order of 1% at the arcminute scale). 

For a reliable comparison of inflationary predictions and CMB data we 
need analytical predictions that can meet the accuracy of the observations. 
The primordial fluctuations are characterised in terms of power spectra of 
scalar and tensor perturbations (inflation driven by a scalar field predicts that 
there are no vector perturbations). Only a limited number of very special 
inflationary models is known for which the power spectra can be calculated 
exactly ||. Thus we have to use some approximation or we have to rely on 
numerical calculations ||. The state of the art of the analytical calculation 
of the spectral amplitudes are the approximate expressions due to Stewart 
and Lyth @ , which are first order expressions in terms of the so-called, slow- 
roll parameters. These expressions have been tested in Refs. H |9j to be 
precise enough to match the accuracy of current observations, given the con- 
ditions for slow roll are met. However, it has been pointed out that for the 
analysis of data from future observations, more precise analytical expressions 
are compulsory ||. Moreover, there are inflationary models that do not be- 
long to the class of slow-roll models. The approach by Stewart and Lyth 
includes several approximations. In a first step the mode equations of 
the fluctuations are approximated by Bessel differential equations. With this 
aim, one has to assume that the slow- roll parameters are constant, which is a 
consistent procedure up to first order in the slow-roll parameters only. A sec- 
ond approximation is made by assuming that the dominant modes take their 
superhorizon values at the moment of horizon crossing, which is arbitrarily 
fixed as k = aH. Martin and Schwarz || closed a gap in the derivation of 
this approximation, showing that the time of matching the Bessel function 
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solution to the superhorizon solution has to coincide with the time when the 
Hubble rate and the slow-roll parameters are evaluated. Wang, Mukhanov 
and Steinhardt JTD[ have shown that the method of Stewart and Lyth [7[ does 
not allow to improve the Bessel function approximation to higher orders. Re- 
cently, Stewart and Gong 0] introduced a method to calculate the spectra 
to higher orders in the slow-roll parameters for slow-roll inflation. This new 
method does no longer refer to the Bessel function approximation, but in- 
stead uses an approximation based on Greens functions. They calculate the 
scalar spectral amplitude at second order in the slow-roll parameters. 

In this paper we present an approximation that does not rely on the slow- 
roll approximation and has the advantage that models for which some of the 
slow-roll parameters take large values can be treated. We calculate the am- 
plitudes with high precision for two families of inflationary models, models 
that are characterised by an almost constant Hubble distance during inflation 
and models that are characterised by an almost linearly growing Hubble dis- 
tance. The common link between both families is that their spectral indices 
are almost constant, even for large values of some of the slow-roll param- 
eters. Our approximation makes use of the Bessel function approximation, 
but avoids the argument by Wang, Mukhanov and Steinhardt |T(J by putting 
constraints on the relative magnitude of the various slow-roll parameters. 

The traditional set of slow- roll parameters (see e.g. [T^]) is not well suited 
for the new approximation. Although we hesitated to introduce yet another 
set of parameters that control the dynamics of inflation, our new parameters 
have two major advantages. Firstly, their definition is simpler than other 
definitions of slow-roll parameters, and, secondly, the new definition allows a 
very transparent physical interpretation of the parameters: the parameters 
control how the Hubble distance behaves during inflation and thus we call 
them the horizon flow functions (or parameters when we evaluate these func- 
tions at a certain moment of time). We have presented a variant of our new 
approach using the traditional notation in Ref. JT3| . 

This paper is organised as follows: In section 2 we introduce the horizon 
flow functions and establish the link to the slow-roll parameters. The core 
of the paper is section 3, where we calculate the scalar power spectrum in 
the constant-horizon and in the growing-horizon approximations up to third 
order. The tensor power spectra are calculated in section 4. In section 5 
we present the corresponding 'consistency relations of inflation'. We finally 
compare our results to the result of Stewart and Gong |TT] and show that 
both results are consistent, which is a nontrivial check of both calculations. 
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2 Horizon flow 



For the reasons explained in the introduction we define a set of horizon flow 
functions starting from 

"Hi 

where du = 1/H denotes the Hubble distance, iV = ln(a/a;) the number 
of e-folds since some initial time ti, and dm = d^(ti). (Note that usually 
the number of e-folds is counted backward in time, we count it forward, i.e., 
N(t[) = 0.) The quantity dn is commonly called "horizon" because it is a 
good estimate of the size of the region that may be in causal contact within 
one expansion time. 

Further we define a hierarchy of functions in a systematic way by 

_ dln|e m | 

Wi = dN , m>0. (2) 
According to this definition, 

dlnd H /0 s 

measures the logarithmic change of the Hubble distance per e-fold of expan- 
sion. Inflation happens for e\ < 1 (equivalent to a > 0) and e\ > from 
the weak energy condition (for a spatially flat universe). For m > 1, e m may 
take any real value. Expressions (0) define a flow in the space {e m } with the 
cosmic time being the evolution parameter. This flow is described by the 
equations of motion 

1 

^o e m " -] — e m e m _t_i — 0. (4) 

"Hi 

For m = 0(1) we find ex = g?h and = c^h^h respectively, which describe 
the time evolution of the horizon. 

Let us stress the advantages of definition (0) over other definitions of the 
slow-roll parameters |7|, || |I^] . First of all, the physical interpretation of the 
functions given by Eq. (Q) and equations of motion (|4]) is straightforward 
and model independent (not restricted to models with a single scalar field). 
Secondly, the notation is concise, leading to significant simplification in the 
involved expressions. Thirdly, the definition is easy to memorise. The link 
between the first three horizon flow functions and various definitions of the 
corresponding slow-roll parameters is presented in table [l]. 
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Lidsey et al. [12 


Martin & Schwarz || 


Stewart & Lyth [7^ 


work 






Stewart & Gong 


ei 
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e 


ei 




2e - 2r] 


2e-2<5 


2ei + 25i 


£263 


4e 2 - 6er] + 2£ 2 


2£ 


4e 2 + 6ei 5i - 25 2 + 25 2 



Table 1: Conversion table for different definitions of the expansion basis. 



a e — > ei and <5 — * r5i 



To proceed with the calculation of the inflationary perturbations, we have 
to relate the comoving Hubble rate, aH, to conformal time. From the defi- 
nition dr = dt/a we find after a partial integration 

1 + /^. (5) 



off(l-ei) ' 7 (1 -ex) 2 aH 

The equation of motion of the perturbations is (see [0, |9[ for the nota- 
tion): 



fl(k,T)" + (k 2 - — )fJ,(k,T) =0, 



(6) 



where a prime denotes a derivative with respect to conformal time and where 
z = dyfei for scalar perturbations and z = a for tensor perturbations. This 
equation should be solved with the following initial conditions 

e -ik(r-Ti) 

lim // s ,t(t) = =f4^pi 7=^, (7) 

k/{aH)^+oo 



where Zpi denotes the Planck length (the two signs stand for scalar and tensor 
perturbations respectively). Then the power spectra can be calculated and 
read 

2k 3 ' 



k 6 Pc 



k 3 




8tt 2 





k 6 P h = — 



(8) 



where ( and h stand for scalar and tensorial modes respectively. 

3 Scalar perturbations 



The potential in the scalar mode equation reads 

Z ry TT r) [ n 3 If) 1 1 

: v 2 



a 2 H 2 ( 2 - ei + ^e 2 + -e\ - -e x e 2 



£263 



(9) 



We can solve the mode equations by Bessel functions 



H = (krY^lBMkr) + B 2 J_ u (kr)] , (10) 
with the constants B\ and B 2 , if 

4^-1 (id 

may be approximated to be constant. This condition is met if we can con- 
sistently neglect the time derivative of e\ and e 2 , which means that we have 
to set e\e 2 ~ and €263 ~ from equations Note that this does not 
imply e n ~ 0. We see from equation (||) that we actually need to require 
that eie 2 AN < A/100%, where A is the accuracy that we want to achieve 
and AN is the number of e-folds during which the time derivatives have to 
be neglected. Dropping the above mentioned terms in Eqs. (|9|) and @ we 
find the index of the Bessel functions, which is given by 

111 

^2 + — + 2 £2 ' (12) 

In the above expression, 1/(1 — ei) should be expanded for small e\ and 
truncated at order n, such that e™ > |eie 2 |. 

We can now fix Bi, B 2 by comparing to Eq. ([?[) in the limit r — > T{. This 
gives 

, exp[z(z/7r/2 + ix/A + krd] , , s 

B 1 = 2nl Pl PU ' -L- ^, B 2 = - exp(-i7ri/)Si. (13) 

V k sm(7T^) 

The next step is to calculate the superhorizon limit of /i, which becomes 

^ / P1 7r2"+ 1 (-fcr) 1 / 2 -" _ l vx 2 v+ \-kT) l / 2 - v Y(u) 
11 Vk sin(7rz/) 

Inserting in the definition of the power spectrum (0) and eliminating r with 
help of Eq. (13) we find 



k 3 P ( ~ P1 



7T 



, v(3-a/) 

^ • (15) 



off 



At this point we have to fix the time t* at which we evaluate e m . We introduce 
the abbreviation fc* = (aH)(t*). We have shown above, that terms like e\e 2 
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should be neglected consistently. However, the above result certainly still 
involves terms like this. In the last step we therefore have to expand Eq. flT5|) 
for e± and e 2 and drop all terms that are of order eie 2 or smaller. The crucial 
difference to the slow-roll approximation as used by Stewart and Lyth |7]] is 
that they assume e = e\ and 5 = €2/2 — t\ to be small simultaneously, thus 
they cannot keep higher power in their slow-roll parameters. In our approach 
we can, e.g., keep terms of order e* m if they are larger than |eie 2 | and |e 2 e 3 |. 
Consequently, depending on the relative magnitude of the functions ei, e 2 
and €3 we obtain different results. 



3.1 Constant-horizon approximation 

In some inflationary models as for inflation at a maximum, the time derivative 
of the Hubble distance is tiny. For this kind of models, during a certain num- 
ber of e-folds, 6\ <C 1. However, it does not necessarily mean that all other 
e m have to be small as well. We thus define the constant-horizon approxima- 
tion at order n for the situation je^l > max(|eie2|, l^^l); which means that 
we are allowed to include the following monomials in the primordial spectra: 
1, ei, e 2 , • • • , e 2 . 



Expanding Eq. flT5|) we find at third order 



k S Pr 



a 



a 1 



a 2 



a + d\ In 



7rei 

l-2(C+l)ei 
4C 3 



+ a 2 In 2 




Ce 2 + -(4C 2 + 7r 2 - 



1 

24 



3C(8 



7T 



+ 14C(3) - 16 



)4 



-2 ! 



-26! 
1 



£2 



Ce 2 2 --(AC 2 + 7i 2 



-2 5 



2 — o^ 6 2 ' 
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a 3 = 

where C = 7e + In 2 — 2 



(16) 

(17) 
(18) 
(19) 
(20) 



-0.7296 and C(3) w 1.2021. The quantity that 
is usually called the amplitude is obtained by setting k = k*. We finally 
calculate the spectral index 



n s - 1 = 



d In k 3 P ( 
d In k 



(21) 
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which is most easily obtained from Eq. It reads n$ — I = — 2e% — e 2 

at any order n in the constant horizon approximation. In the corresponding 
limits this result agrees with the usual slow-roll expression. Consistently with 
our assumptions, there is no "running" of the spectral index. 



3.2 Growing-horizon approximation 

Power-law inflation (a oc t p ) is one of the few models of inflation for which an 
exact expression for the primordial spectrum can be obtained in closed form. 
Our horizon flow functions in this case are ei = 1/p = const, and e m = for 
m > 1. This observation suggests that inflationary models with \e m \ < €\ for 
m > 1 might be approximated by keeping all terms in e\ up to order n, where 
n is the maximal integer for which e™ > max(|eie2|, l^sl) holds true. For 
this kind of models, ei ~ const., implies that the Hubble radius grows almost 
like a linear function of time. Hence, we define the (linearly) growing-horizon 
approximation to order n to include the following terms: 1, ei, . . . , e™, e 2 . 

The scalar power spectrum at third order in the growing horizon approx- 
imation reads: 



k 3 P, 



b 



bi 



7rei 



61 In (^-)+h\n 2 (^-)+h\n* (f ] + . . . 



k 



k 



1 



1 - 2(C + l)ei + - [AC{C + 1) + 7T 2 - 10 



1 r 
3 



AC 3 + 3C(tt 2 - 12) + 14C(3) - 19 e\ - Ce 2 

-2 Cl + 2(2C + l)t\ - {AC 2 + 7T 2 - \2)e\ - e 2 , 
2e 2 - ACe\ , 
4 , 



— e 



i- 



(22) 



(23) 

(24) 
(25) 

(26) 



Ignoring the monomials e 2 and ef, but keeping e 2 , we recover the result of 
Stewart and Lyth M. The spectral index at n-th order reads 



ns - 1 



2u 



-2(e 1 + e 2 1 + el + ... + e^)-e 2 , 



(27) 



and again there is no "running" of the spectral index. 
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4 Tensor perturbations 



Now the potential in the mode equation reads 

(aHf (2 - ei ) , 



z 
z 



(28) 



which allows us to solve the mode equation with the same approximation as 
in the scalar case, but with Bessel function index 



v 



1 
2 



(29) 



The following steps are analogous to those for the calculation of the scalar 
perturbations, the basic difference is that no contribution of e 2 arises in any 
case. 

Using the constant-horizon approximation we obtain 



k 3 P h 



a 
ai 



16l 2 Pl H 2 



TV 

1-2(C 

-2e x , 

-2e x , 



a ° + ailn + 
- Dei, 



(30) 

(31) 
(32) 
(33) 



at any order. Note that all aj for i > 2 vanish in the constant horizon 
approximation and that there are no higher order corrections to n^, which 
follows from the absence of any e2-dependence in the potential ( |28"D and the 
fact that higher powers of €\ are consistently neglected. The growing-horizon 
approximation leads to 



k 3 P h 



m 2 m H 2 

7i 



b + bi In 



1 



bo In 2 



+ b 3 In 3 



1 - 2(C + l)ei + - [4C(C + 1) 



tt- 
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1 



AC 3 + 3C(vr 2 - 12) + 14((3) - 19 



bi 
b 2 

b 3 



-2e a + 2(2C + l)e\ - {AC 2 + tt 2 - 12)e 3 , 

2e\ - 4Cel 
4 



-3* 



(34) 



(35) 

(36) 
(37) 

(38) 



9 



at third order, and 



nT = _ 2 ( ei + e 2 + e 3 + ... + <£), (39) 

at nth order. Again both results are consistent with the results of Stewart 
and Lyth 0. 



5 Consistency relations 

It is interesting to inspect the so-called 'consistency relations of inflation' 
for models that do not belong to the class of slow-roll inflation. We define 
the tensor-to-scalar ratio r = Ph/P(. The 'classic' result for slow-roll models 
reads at first order tit = — r/8. We find for the constant horizon and for the 
growing horizon approximations r = 16ei at any order. This result may be 
used to express the corresponding expressions for the tensorial spectral index 
as 

n T = - T - (40) 



for the constant horizon approximation at any order, and 



(41) 



for the growing horizon approximation at order n. Both results are consistent 
with the second order result of Ref. |TJ], = — 2[(r/16) — (r/16) 2 — (ns — 
l)(r/16)], since (ns — l)(r/16) = — 2(r/16) 2 + C(eie 2 ). For models where 
the constant horizon approximation applies we expect that r < 1, thus the 
prospects to detect the tensors and therefore to test the consistency relation 
are bad. In the case of those models for which the growing horizon approx- 
imation is suited, there is a chance to detect the tensor contribution. Here, 
the corrections to the slow-roll results given by Eq. ( f4T|) may be relevant. 
However, current data seem to indicate that r < 1, which implies that higher 
order corrections to the classic consistency relation are not important. 

6 Conclusion 

We may compare our result with the work by Stewart and Gong [[□]]. They 
used a new method to obtain the scalar amplitude at second order in the 
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Figure 1: Regions in the ei- [ 62 1 parameter space where the spectral am- 
plitudes can be calculated with an accuracy better than 1%. In the dark 
shaded region the Stewart-Lyth (SL) approximation as well as all other 
approximations are fine. Second-order corrections, as calculated by Stewart 
and Gong (SG) JTl|], extend that region to the light shaded region. The 
constant horizon approximation at order n (elm), and the growing horizon 
approximation at order n (ghn), do well below the thick line. The rays in- 
dicate where the corresponding higher order corrections are necessary The 
thick line itself is the condition e 1 |e 2 | < (A/ 100%) /AN, with AN = 10 and 
A = 1%. As is easily seen the ch2 and gh2 regions are included within the 
SG region. The chn and ghn regions with n > 2 allow us to go beyond the 
SG approximation. 
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slow-roll parameters. In our notation they obtain: 



k 3 P, 




e 2 e 3 ,(42) 



which agrees with our result if eie 2 ~ and e 2 e 3 m 0. 

Our results can be applied to many interesting models of inflation. The 
constant horizon approximation works fine if the inflaton field sits close to a 
maximum of the potential. On the other hand the growing horizon approxi- 
mation is a good approximation for chaotic inflation models with monomial 
potential oc <p p , when p > 4. To give an example, for p=8we have e\ ~ 0.04 
and e 2 ~ 0.02, which allows us to take all quadratic terms in t\ into account. 
We plot in figure 1 the domain of applicability of the various approximations 
in the ei- 1 62 1 plane. We tacitly assume that \e^\ < min[ei, | €2 1 ] , for the sake 
of simplicity of the argument. If this condition is not satisfied an analogous 
three dimensional plot has to replace figure 1. 
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